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CHANDRA OBSERVATION OF PSR B 1823- 13 AND ITS PULSAR WIND NEBULA 
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ABSTRACT 

We report on an observation of the Vela-like pulsar B 1823- 13 and its synchrotron nebula with the Chandra 
X-ray Observatory. The pulsar's spectrum fits a power-law model with a photon index TpsR ~ 2.4 for the 
plausible hydrogen column density «h = 1 X 10 22 cm" 2 , corresponding to the luminosity Lpsr « 8 x 10 31 ergs 
s" 1 in the 0.5-8 keV band, at a distance of 4 kpc. The pulsar radiation likely includes magnetospheric and 
thermal components, but they cannot be reliably separated because of the small number of counts detected 
and strong interstellar absorption. The pulsar is surrounded by a compact, ~ 25" x 10", pulsar wind nebula 
(PWN) elongated in the east- west direction, which includes a brighter inner component as 7" x 3", elongated 
in the northeast-southwest direction. The slope of the compact PWN spectrum is r comp w 1.3, and the 0.5-8 
keV luminosity is L comp ~ 3 x 10 32 ergs s" 1 . The compact PWN is surrounded by asymmetric diffuse emission 
(extended PWN) seen up to at least 2.4' south of the pulsar, with a softer spectrum (r ext w 1.9 for «h = 1 x 
10 22 cm" 2 ), and the 0.5-8 keV luminosity L ext ~ 10 33 -10 34 ergs s . We also measured the pulsar's proper 
motion using archival VLA data: [i a = 23.0 ±2.5 mas yr" 1 , fig = -3.9 ±3.3 mas yr" 1 , which corresponds 
to the transverse velocity vj_ ~ 440 km s" 1 . The direction of the proper motion is approximately parallel to 
the elongation of the compact PWN, but it is nearly perpendicular to that of the extended PWN and to the 
direction towards the center of the bright VHE 7-ray source HESS J 1825-1 37, which is likely powered by PSR 
B1823-13. 

Subject headings: ISM: individual (G18.0-0.7, HESS J1825-137) — pulsars: individual (PSR B1823-13 = 
J 1826-1 334) — stars: neutron — X-rays: ISM 



1, INTRODUCTION 

About a dozen of ~ 80 radio pulsars observed in X-rays be- 
long to the group of the so-called Vela-like pulsars, with spin- 
down ages r = P/2P = 10-30 kyr and spindown powers E of 
a few xlO 36 ergs s" 1 . Their radiation usually show a ther- 
mal component, possibly emitted from the cooling neutron 
star (NS) surface, and a nonthermal component with a power- 
law (PL) X-ray spectrum, presumably originating in the pul- 
sar magnetosphere (see Kaspi et al. 2006 for a review). These 
pulsars are enveloped by pulsar wind nebulae (PWNe) whose 
X-ray emission is due to the wind of ultrarelativistic particles 
shocked in the ambient medium (Kargaltsev et al. 2007, and 
references therein). Therefore, in addition to studying pul- 
sar magnetospheres and NS cooling, X-ray observations of 
Vela-like pulsars allow one to study pulsar winds, the ambi- 
ent medium, and their interaction. 

PSR B 1823-1 3 (also known as PSR J 1826-1 334, hereafter 
referred to as B1823) is a typical Vela-like pulsar. It was de- 
tected in the Jodrell Bank radio survey (Clifton et al. 1992). 
Its period, P = 104 ms, and period derivative, P = 7.5 x 10~ 14 s 
s , imply the spindown age r = 21 kyr, spindown power E = 
2.8 x 10 36 ergs s" 1 , and surface magnetic field = 2.8 x 10 12 
G. The distance to the pulsar, d rj 4 kpc, was estimated from 
the pulsar's dispersion measure (DM = 231 pc cm" 3 ) using 
the models of Galactic electron distribution (Taylor & Cordes 
1993; Cordes & Lazio 2002). Although such a young pulsar 
is expected to reside in a supernova remnant (SNR), observa- 
tions of this field with the Very Large Array (VLA) have not 
found an SNR, nor a radio PWN (Braun, Goss, & Lyne 1989; 
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Frail & Scharringhausen 1997; Gaensler et al. 2000). 

This pulsar is particularly interesting because it possibly 
powers an extended (~ 1°, i.e. ~ 70 pc, in size) Very High 
Energy (VHE) 7-ray source HESS J 1825-1 37 (Aharonian et 
al. 2005,2006). The 7-ray emission extends asymmetrically 
south-southwest of B1823, with the peak surface brightness 
w 10' from the pulsar. Aharonian et al. (2006, hereafter A06) 
found a softening of the 0.25-35 TeV spectrum with increas- 
ing distance from the pulsar (hence an energy-dependent mor- 
phology), with a photon index T 7 increasing from ks 1 .9 to 
« 2.5. These authors discuss a scenario in which the 7- 
rays are produced by Compton upscattering of the cosmic 
microwave background (CMB) photons off the ultrarelativis- 
tic electrons supplied by the pulsar. As the 7-ray luminos- 
ity of the source, L 7 ~ 3 x 10 35 af| ergs s" 1 above 200 GeV 
[c/4 = of/(4kpc)], is only a factor of 10 lower than the pulsar's 
spindown power, A06 propose that the spindown power was 
significantly higher in the past, when the electrons responsi- 
ble for the currently observed TeV emission were produced 
by the pulsar. 

Although the field of B1823 has been observed with 
ROSAT, ASCA, and XMM-Newton, the X-ray spectrum of 
the pulsar and the structure of its surroundings remained un- 
clear because of the relatively low angular resolution of those 
missions. Finley et al. (1996, F96 hereafter) observed B1823 
with the ROSAT observatory and, despite the small number 
of source counts detected (e.g., 127 counts from a 5' radius 
circle around the pulsar position in the 14.7 ks PSPC expo- 
sure) found evidence of a point source surrounded by a com- 
pact PWN of ~ 20" radius, and a dimmer component ex- 
tending up to ~ 4' south-southwest of the pulsar. Because 
of the poor statistics, it was impossible to measure the spec- 
tra. For an upper limit on the hydrogen column density, 
«h.22 = nu/ (10 cm" 2 ) < 4, F96 estimated the upper limit 
fluxes at a few xlO" 12 ergs cm" 2 s" 1 level for the emission 
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components, corresponding to luminosities < (4-8) x 10 33 <f| 
erg s -1 in the 0.5-2.4 keV range. Sakurai et al. (2001) were 
able to detect the extended structure around the B1823 posi- 
tion with ASCA, but that observation lacked the spatial reso- 
lution to separate the components. 

Gaensler et al. (2003, G03 hereafter) observed B 1823 with 
the XMM-Newton EPIC detectors. The XMM-Newton reso- 
lution turned out to be too low to resolve the pulsar, but G03 
detected two components of the PWN (designated as G18.0- 
0.7 based on its Galactic coordinates) associated with this 
pulsar: a "core" of a 30" extent, elongated in the east-west 
direction, and a large-scale diffuse component of lower sur- 
face brightness, ~ 5' in extent, seen mostly to the south of the 
pulsar. Fitting the core and diffuse component spectra with 
a PL model, G03 found the photon indices r core « 1.6 and 
r dif f w 2.3, and 0.5-10 keV luminosities L core w9x I0 32 dj 
ergs s" 1 , Ldiff ~ 3 x 10 33 <f| ergs s" 1 , for a foreground hy- 
drogen column density n H ,22 = 1 -2. They suggested that the 
core represents the pulsar wind termination shock at a dis- 
tance R s w 15" from the the pulsar while the diffuse com- 
ponent indicates the shocked downstream wind. G03 pro- 
pose that the asymmetric morphology of the diffuse com- 
ponent with respect to the pulsar is the result of a reverse 
shock from an associated SNR, which has compressed and 
distorted the PWN. They place limits on the pulsar's ther- 
mal emission: kT < 0.147 keV and bolometric luminosity 
Lboi < 8.7 x 10 33 c/| ergs s" 1 (as seen by a distant observer), 
assuming that the NS emits as a blackbody of a radius R = 12 
km, for « H ,22 < 2. 

Because of its superior angular resolution, Chandra is suited 
much better than any other X-ray observatory for studying 
small-scale extended structures (such as the core of the B 1823 
PWN) and resolving point sources possibly projected onto 
such structures. Therefore, our group proposed a Chandra 
observation of B1823, with the goal to separate the pulsar 
from the PWN, measure their spectra, investigate the PWN 
morphology with the best possible resolution, and compare 
the pulsar and PWN properties with those of other Vela-like 
pulsars. The data from that observation immediately showed 
that the pulsar can indeed be separated from the PWN, the 
"core" PWN is highly nonuniform in surface brightness, and 
the overall PWN morphology is more complicated than it had 
been presumed based on the previous low-resolution X-ray 
observations. These results have been reported at several con- 
ferences (Teter et al. 2003a,b; Gonzalez & Kaspi 2006; Kar- 
galtsev et al. 2006), but they have not been published in a 
refereed journal. In this paper, we decribe the analysis and 
results of the Chandra observations and present infrared (IR) 
and radio images of the field (§2). As the interpretation of the 
PWN morphology strongly depends on the direction and mag- 
nitude of the pulsar's proper motion, we also report the proper 
motion measurement in §3. Implications of our findings for 
the pulsar and PWN physics, including the VHE 7-radiation, 
are discussed in §4, while §5 summarizes the results of this 
work. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. Chandra ACIS observation 

B1823 was observed with the Advanced CCD Imaging 
Spectrometer (ACIS) aboard Chandra on 2002 October 24 
(ObsID 2830). Although last 13 ks of the 40,836 s exposure 
suffered from flaring particle background (9.7 x 10~ 7 counts 
s" 1 per pixel, on average), we have checked that filtering out 



the periods of elevated background at the expense of short- 
ening the exposure time would actually reduce the signal-to- 
noise ratio for the pulsar and the compact PWN and only very 
slightly increase it for the large-scale extended PWN; there- 
fore, we chose to use the whole exposure for our analysis. 

The observation was carried out in Faint mode, and the tar- 
get was imaged on S3 chip, with the radio pulsar position 
« 35" from the optical axis, where the angular resolution is 
still not degraded. The other ACIS chips activated during this 
observation were SI, S2, S4, 12 , and 13. The detector was 
operated in Full Frame mode, which provides time resolution 
of 3.24 seconds. We analyzed the data reprocessed by the 
Chandra X-ray Center (CXC) on 2006 October 1 1 (rev. 2; AS- 
CDSVER 7.6.9). The data were reduced using the Chandra 
Interactive Analysis of Observations (CIAO) software (ver. 
3.4; CALDB ver. 3.3.0.1). We used the software package 
XSPEC (ver. 1 1.3.2) for the spectral analysis. 

2.2. X-ray images ofB1823 and its vicinity 

Figure 1 shows the ACIS-S3 images of the B1823 field, 
with an extended source around the radio pulsar position. 
The coordinates of the highest peak in the X-ray bright- 
ness distribution (which is at the center of a 2 x 2 pixel is- 
land, with equal numbers of counts in each of the 4 pix- 
els) are a= 18 h 26 m 13. 187 s , <5 = -13°34'46.63" (J2000), with 
an estimated la centroiding error of 0.2" for each of the 
coordinates. This position differs by 0.08" in right ascen- 
sion and 0.17" in declination from the radio pulsar position, 
a = 18 h 26 s 13.182 s ±0.04", S = -13°34'46.8" ±0.2", calcu- 
lated from the original position in Hobbs et al. (2004) for the 
epoch of the Chandra observation (MJD 52,571), using the 
proper motion measurement described in §3. These differ- 
ences are smaller than the combined errors (0.3" in right as- 
cension and 0.4" in declination, at the 68% level) of the peak 
position measurement, radio position, and Chandra absolute 
astrometry. Therefore, there are no doubts that the highest 
brightness peak corresponds to the pulsar B1823 while the 
surrounding emission represents its PWN. 

We see from Figures 1 and 2 that the PWN morphology 
is asymmetric and nonuniform. The brightest inner PWN 
component (surface brightness in the range of 4-23 counts 
arcsec -2 in 0.7-7 keV), is - 7" x 3" (i.e. 0.14 x 0.06 pc 2 at 
d = 4 kpc) in size. It is approximately centered on the pulsar 
and elongated along the northeast-southwest direction (posi- 
tion angle P.A. ~ 50°, measured east of north). 

The inner component is embedded in a dimmer (average 
surface brightness « 1.2 counts arcsec" 2 in 0.7-7 keV) outer 
component of the compact PWN, within a 25" x 10" area. It 
is stretched along the east-west direction, with the pulsar be- 
ing closer to its east end (i.e. the outer component is more 
extended in the direction approximately opposite to the di- 
rection of the pulsar's proper motion; see §3). Its surface 
brightness distribution is rather nonuniform, with a promi- 
nent enhancement w 8" west of the pulsar (we designate it 
as the "blob" and mark 'b' in the figures). The smoothed 
subpixel-resolution image shown in Figure 1 (top right) sug- 
gests that the blob is connected to the inner PWN component 
by a faint "bridge". We have checked that there are no optical- 
NIR counterparts to the blob in the catalogs available (e.g., the 
nearest 2MASS 3 counterpart is 6.5" away). 

Another compact source is seen in the images w 18" east 
of the pulsar. This object, CXOU 182614.4-133448 (here- 

3 Two Micron All Sky Survey (Skrutzkie et al. 2006). 
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FIG. 1.— To/) /e/f: 77" X 48" ACIS-S3 image of B1823 and its PWN (0.7-7 keV; pixel size 0"49). Top right: Adaptively smoothed image (0.7-7 keV; pixel 
size 0"49) of the same region obtained by removing the pipeline pixel-randomization and applying the sub-pixel resolution tool (Tsunemi et al. 2001; Mori et al. 
2001). Bottom left: The same image as in the top left panel showing the extraction regions used for the spectral analysis of the PWN and the pulsar (see §2). 'b' 
and 'X' mark the "blob" in the compact PWN and a point source unrelated to the PWN, respectively. Bottom right: Heavily binned (pixels size 3"94; 0.5 — 8.0 
keV) and smoothed (with a gaussian kernel of r — 27") ACIS-S3 image of B 1823 and its surroundings. The brightness and smoothing scales are chosen to stress 
the fainter, more extended emission south of the compact PWN. The polygon shows the region from which we extract counts for the spectral analysis of the 
large-scale extended emission, with the background estimated from the r — 1 .28' circle northwest of the PWN. The annulus around the compact PWN is used for 
measuring the compact PWN background. 



after referred to as Source X) is point-like. The analysis of 
its spectrum (§2.5) shows that it is likely a background AGN, 
unrelated to B 1823. 

The smoothed image in Figure 2 (top) shows marginal ev- 
idence of two ~ 20"-long "tails", apparently emanating from 
the pulsar in the west-northwest and west-southwest direc- 
tions. The linear profiles in the middle and bottom panels 
of Figure 2 show that the outer PWN structure is discernible 
above the local background level of ~ 0.2 counts arcsec -2 up 
to w 7" and w 23" east and west of the pulsar, respectively. 

The compact PWN (comprised of the above-described 
inner and outer components) is, in turn, immersed in an 
even dimmer large-scale diffuse emission apparent in heav- 
ily binned images. As it is natural to assume that this diffuse 
emission is also produced by the wind from B1823, we will 
call it the extended PWN. It is mostly concentrated southward 
of the compact PWN, and it is discernible in the ACIS im- 
ages up to the edge of the S3 chip (i.e. « 2.4' from the pulsar; 
Fig. 1, bottom right). (It likely extends further into the S2 
chip, but the lower sensitivity of this chip does not allow a re- 
liable analysis of the faint diffuse emission.) The extended 
PWN is somewhat better seen in the XMM-Newton EPIC 
MOS1+MOS2 image shown in Figure 3c (see also G03). 
However, the true extent of this dim diffuse component is dif- 
ficult to determine even in the large field-of-view of the MOS 
detectors because of the strong vignetting at large off-axis an- 



gles and high EPIC background. 

2.3. PWN spectrum 

We extracted the compact PWN spectrum 4 (Fig. 4, top) 
from the ellipse shown in Figure 1 (bottom left), excluding 
the counts within the r = 3" circle centered on the X-ray po- 
sition of the pulsar (see §2.1) and restricting the energy range 
to 0.7-7 keV. This region of 225 arcsec 2 area contains 415 
counts. We estimated the background from the annulus of 
2916 arcsec 2 area shown in Fig. 1 (bottom right), which con- 
tains 511 counts. Thus, there are 375 ±20 compact PWN 
counts in the aperture chosen. 

To analyze the spectrum of the extended PWN (Fig. 4, 
bottom), we chose the polygon encompassing the region of 
brightest emission south of the pulsar, with the area of 3.4 
arcmin 2 (see Fig. 1, bottom right), which contains 2329 
counts in the 0.7-7 keV band. The background was extracted 
from the r = 1.28' circular region shown in the same panel, 
which contains 2004 counts. Subtracting the background con- 
tribution, we obtain 1008 ±57 extended PWN counts in the 
polygon aperture. (We have also tried several other back- 
ground regions but this did not lead to significant changes in 
the spectral fits.) 

4 As we found no statistically significant differences between the spectra 
of the inner and outer components of the compact PWN, we present here only 
the spectrum of the entire compact PWN. 
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FIG. 2.— Top: 62" X 48" ACIS-S3 image (0.7-7 keV; pixel size 
0.49"; smoothed with 1.5" gaussian kernel) of the immediate surroundings 
of B1823, including the compact PWN and Source X (likely a backlground 
AGN). The white arrow shows the direction of the pulsar's proper motion 
(see §3). The length of the arrow, 7", corresponds to the distance that the 
pulsar travels during 300 years. The middle and bottom panels show the 
linear profiles extracted by moving the rectangular box (0.49" X 17" and 
0.98" X 17.1", respectively) along the horizontal line shown in the top panel. 
The one-dimensional point spread function centered on the pulsar position is 
shown in the middle panel. 



First, we fit the compact and extended PWN spectra with 
the absorbed power-law (PL) model, allowing the hydrogen 
column density to vary. The confidence contours obtained 
from these fits are shown in Figure 5. We see that the 90% 
confidence contours for the compact and extended PWN do 
not overlap, which suggests that the two spectra are dif- 
ferent. Since we do not expect the hydrogen column den- 
sity, «h,22 ~ 1, to differ significantly on such a small, ~ 2', 
scale, we ascribe the difference to different spectral slopes, 
Text r comD ~ 0.5. 



To estimate the spectral slopes more accurately, it would be 
useful to restrict the range of «h (which would also be help- 
ful for the spectral analysis of the pulsar - see §2.4). The fits 
with free «h give «h,22 = 1 -24 ± 0.26 and «h.22 = 0.72 ± 0.15 
for the compact and extended PWNe, respectively (the uncer- 
tainties correspond to the 68% confidence level for a single 
interesting parameter). The values of «h.22 obtained by G03 
using the XMM-Newton data are in the range of 0.8— 1.7, de- 
pending upon the spectral model and extraction region cho- 
sen. The estimate based on the pulsar's dispersion measure 
(DM= 231 cm" 3 ) and the usual assumption of 10% degree of 
ionization of the interstellar medium (ISM) gives «h.22 = 0.7. 
This is close to the lower limit, «h,22 ^ 0.7, estimated from the 
observed colors of the foreground (d = 2.9 kpc) 06 star HD 
169727, projected w 8' southeast of the pulsar (F96). Finally, 
the total Galactic HI column in that direction is «hi,22 = 1 -4- 
1.5 (Dickey & Lockman 1990). Since the «h values deduced 
from an X-ray spectrum under the assumption of standard el- 
ement abundances usually exceed the «hi measured from 21 
cm observations by a factor of 1.5-3 (e.g., Baumgartner & 
Mushotzky 2005), this «hi value does not contradict to those 
estimated from the X-ray fits and the dispersion measure. 
Thus, with allowance for different kinds of uncertainties, it 
seems reasonable to assume 0.7 < «h.22 ^ 1-3, with «h.22 ~ 1 
being the most plausible estimate. 

With the hydrogen column density fixed at «h.22 = 1 0, the 
absorbed PL fits to the compact and extended PWN yield pho- 
ton indices r comp ~ 1.3 and r ext ~ 1.9, and 0.5-8 keV lumi- 
nosities L comp ~3 x 10 32 and L ext ~9 x 10 32 ergs s" 1 (see also 
Fig. 6 and Table 1 ; the extended PWN luminosity is for the 
polygon region). 

2.4. Pulsar spectrum 

Although the peak of the surface brightness distribution is 
clearly dominated by the pulsar, the surface brightness of the 
inner PWN grows toward the pulsar's position, which intro- 
duces some ambiguity in separating the pulsar contrubution 
from that of the PWN and complicates the measurement of 
the pulsar spectrum. However, based on the surface bright- 
ness distribution near the peak (see, e.g., the one-dimensional 
profiles in Fig. 2), it seems reasonable to assume that the PWN 
surface brightness reaches a plateau of about 3" in size, cen- 
tered on the pulsar, which can be used to estimate the local 
background at the pulsar position. Therefore, we extracted the 
pulsar spectrum from a small circular aperture (small black 
circle in Fig. 1, bottom left) with the radius of 1 ACIS pixel 
(~ 0.49", 55% encircled energy radius) and took the back- 
ground from the 0.74" < r < 1.5" circular annulus centered 
on the pulsar (shown by the white circles in Fig. 1, bottom 
left). The background contributes 22 counts to 62 counts 
extracted from the source region. Given the small number of 
source counts and the appreciable background contribution, 
we chose not to subtract the background channel by chan- 
nel but rather fit the background spectrum (modeled as an 
absorbed PL model with the same «h as for the source) si- 
multaneously with the source spectrum. 

To obtain better constrained fits with the small number of 
counts available, we are forced to freeze the hydrogen col- 
umn density. Since the interpretation of the pulsar spectrum 
(e.g., thermal vs. nonthermal) may depend on the choice of 
the « H value, we fit the pulsar spectrum with the absorbed PL 
model for three plausible values of the hydrogen column den- 
sity: «h,22 = 0.7, 1.0, and 1.3 (see Table 2 and Fig. 7). As 
expected, the inferred photon indices depend on the choice 



Chandra Observation of PSR B 1823-1 3 



5 



MIPS 24 micron 






IRAC 8 micron 



FIG. 3. — Large-scale images of the B1823 vicinity and the central part of HESS J 1825 observed with different telescopes and at different wavelengths. All 
eight panels show the same area on the sky. The TeV brightness contours (adopted from A06) are overlayed on top of the images in panels b—f. a: AOS S3, S2, 
SI, 12 and 13 image (0.7—7 keV; pixel size 3"94), corrected for the exposure map nonuniformity. b: The same image with point sources removed (except for 
B1823) and gaussian (15.7" kernel) smoothing applied, a XMM-Newton EPIC MOS1+MOS2 combined image (0.5-8 keV; pixel size 10"; smoothed with the 
40" gaussian kernel), d: Spitzer IRAC 8 ^m image from the GLIMPSE survey, e: Spitzer MIPS 24 /xm image from the GLIMPSE survey. /■ NRAO VLA Sky 
Survey (NVSS) image at 1.4 GHz. The star in panels b and c marks the location of the brightest part of the extended IR and radio source A seen in the Spitzer 
and NVSS images. The diamond in panels d—f marks the B1823 position. The small white circle marks the position of the brightness peak in the distribution 
of the extended TeV emission from HESS J1825. The projected place of birth of the pulsar estimated from the measured proper motion (see §3) and pulsar's 
characteristic age (r = 2 1 kyr) is shown by the dashed circle. 

verges to a reasonable (although still poorly constrained) val- 
ues of the surface temperature kT ~ 100 eV, projected emit- 
ting area A ~ 500rf| km 2 (corresponding to the bolometric lu- 
minosity Lboi = 4A<tT 4 ~2 x 10 33 c/| ergs s" 1 ), and the slope 
of the PL component Tpsr ~ 1.9. The fitting parameters be- 
come better constrained if we fix the area, but the best-fit tem- 
perature and luminosity depend on the value of A because of 
the strong A-T correlation. For instance, kT increases from 
97 to 139 eV, and L b p f R decreases from 1 .8 x 10 33 to 3.0 x 10 32 
ergs s" 1 , when A decreases from 500 km 2 (approximate pro- 
jected area of the NS surface) to 20 km 2 (projected area in the 
PL+BB fit for the Vela pulsar spectrum). Given the errors of 
these fits and the uncertainty of ISM absorption, the obtained 
BB parameters are just crude estimates. We can, however, 
treat Lboi ~2 x 10 33 c/| ergs s" 1 as an upper limit for the bolo- 



of n H , changing from r P s R = 2.0 ±0.4 to 3.0 ±0.6 with in- 
creasing «h (see Table 2 for details). The pulsar's absorbed 
flux is F PSR = (1.6 ±0.4) x 10" 14 ergs cm" 2 s" 1 in the 0.7-7 
keV band (corrected for the background and finite aperture). 
The unabsorbed pulsar luminosity in the 0.5-8 keV band in- 
creases from Lpsr « 0.6 x 10 32 to 1.1 x 10 32 ergs cm" 2 s" 1 , 
respectively. 

The steepness of the best-fit PL spectra, especially in the 
case of larger tin, indicates the presence of thermal emission 
from the NS surface, such as seen in other Vela-like pul- 
sars (e.g., Pavlov et al. 2001). However, adding the second 
blackbody (BB) component is not justified statistically (i.e. 
the contribution of the BB component is negligible, and the 
fitting parameters of the BB component are unconstrained) at 
nil 22 = 0.7 and 1 .0, while in the case of «h 22 = 1 -3 the fit con- 
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FIG. 4. — Spectrum of the compact PWN (top) and the extended PWN 
(bottom) fitted with the PL model at fixed «h.22 = 1 -0. 
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FIG. 5. — Confidence contours (68% and 90%) in the rtfi-T plane for the 
PL fit to the compact and extended PWN spectra. The contours are obtained 
with the PL normalization fitted at each point of the grid. 

metric luminosity. 

2.5. Other X-ray sources 

Among X-ray sources serendipitously detected in this ACIS 
observation, most interestng are CXOU 182614.4-133448 
(Source X, see Figs. 1 and 2) w 18" east of the B1823 pulsar, 
and a few sources on the S2 and 12 chips projected within the 
brightest part of HESS J 1825 (see Fig. 8), which might supply 
additional ultrarelativistic electrons to power HESSJ1825. 
The positions and spectral characteristics of these sources are 
presented in Table 3. 

The strongly absorbed («h.22 ~ 3) point-like Source X 
shows a hard spectrum (T ~ 0.2) and the 1-8 keV flux of 
~ 3 x 10~ 14 ergs cm" 2 s" 1 (see Table 3 for other proper- 
ties). There is no an obvious optical-NIR counterpart for this 
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FIG. 6. — Confidence contours (68% and 90%) for the PL fit to the compact 
PWN and extended PWN for fixed «h.22 = 1-0. The PL normalization is in 
units of 10~ 5 photons cnT 2 s keV~' at 1 keV. The dashed curves are the 
loci of constant unabsorbed flux in the 0.5-8 keV band; the flux values near 
the curves are in units of 10~ 13 " r "° o-l 
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FIG . 7 . — Confidence contours (68% and 90%) for the PL fit to the pulsars's 
spectrum. The PL normalization is in units of 10~ 5 photons cm -2 s keV~' 
at 1 keV. The dashed curves are the lines of constant unabsorbed flux in the 
0.5-8 keV band (the flux values are in units of 10~ 14 ergs cirT 2 s -1 ). 



source in stellar catalogs. 2MASS images show a point source 
with J = 13.4, H = 12.3 and K„ =11.85 close to Source X 
(a = 18 h 26 m 14.32 s , 6 = -13°34'48.2"), but, since the offset 
of 1 .6" from the pulsar position exceeds the combined cen- 
troiding and absolute astrometry error (0.7" at the 90% confi- 
dence level), the 2MASS source is unlikely to be a counterpart 
of Source X. Based on the X-ray spectrum, we conclude that 
Source X is, most likely, a background AGN, not related to 
either B 1 823 or HESS J 1 825 . 

The other X-ray sources listed in Table 3 also look point- 
like (the apparent extension of the sources s4 and il can be 
explained by the asymmetric PSF broadening at the large off- 
axis distances). 2MASS and Digital Sky Survey (DSS 5 ) im- 
ages do not show NIR-optical objects within most conser- 
vative error circles around the sources si, s5, and il (the 
nearest objects are offset by 4.1", 4.0", and 3.4", respec- 
tively). The strongly absorbed si and s5 are likely extra- 
galactic sources (possibly AGNs), while the nature of il is 

5 see http://archive.eso.org/dss/dss 
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FIG. 8. — ACIS field sources in the vicinity of the HESS J1825 center 
(shown by the circle). 



unclear. The source s4 shows a possible NIR-optical coun- 
terpart DENIS J 182620.9- 134425 (offsets are in the range of 
0.36"-0.86" in different catalogs; magnitudes K s = 11.75, 
H = 12.2, J = 13.3 in the 2MASS catalog; / = 15.0, R = 
15.8-16.4, V = 16.5-17.1, B = 17.5-19.1 in other catalogs). 
The source looks slightly extended in the NIR-optical im- 
ages, likely a blend of two field stars. Close to the sources 
s2 and s3, there are point-like NIR-optical objects (DENIS 
J 1 826 1 7. 1- 1 34 1 1 1 and J 1 826 1 8 .0- 1 34225, respectively), but 
their relatively large offsets from the X-ray positions (1.4"- 
1.9" and 1 .8"— 2.3", respectively, somewhat different in dif- 
ferent catalogs) make their association with the X-ray sources 
questionable. Overall, the X-ray sources in the brightest part 
of HESS J 1825 are rather faint and unremarkable, and there 
is no evidence that any of the sources is related to the VHE 
7-ray source HESS J1825. 

2.6. Vicinity of PSR B1823-13 and HESS Jl 825-1 37 at other 
wavelengths 

To check whether or not there are luminous sources within 
the central part of HESS J 1825, which could enhance the 
VHE 7-radiation by providing additional photons for the in- 
verse Compton (IC) scattering off ultrarelativistic electrons, 
we have examined optical-NIR, IR, and radio images of this 
field. Three examples, Spitzer IRAC 8 /im and MIPS 24 /im 
images from the Galactic Legacy Infrared Mid-Plane Survey 
(GLIMPSE 6 ) and a NRAO VLA Sky Survey (NVSS; Con- 
don et al. 1998) 1.4 GHz image are shown in Figure 3. The 
images do not show anything remarkable in the immediate 
vicinity, within ~ 2', of B1823. The most interesting object 
within the brightest part of HESS J1825 is the extended (« 4' 
in size) source, marked A in Figure 3, prominent in both the 
IR and radio images. The source with a two-sided shell mor- 
phology is centered at a w 18 h 26 m 02 s , 5 w -13°38'08" (mea- 
sured in the IRAC 8 image), i.e. between the pulsar and 
the brightest part of HESS J1825. It is listed in the IRAS cat- 
alog of point sources 7 as IRAS 18231-1340, with fluxes of 
4.25, 13.4, 161, and 393 Jy at 12, 25, 60, and 100 /um, respec- 
tively. Its radio fluxes are 0.13 and 0.03 Jy at 4.85 and 1.4 
GHz, respectively (Griffith et al. 1994; Condon et al. 1998). 

6 http://www.astro.wisc.edu/sirtf/ 

7 Infrared Astronomical Satellite Catalog, 1988, The Point Source Catalog, 
ver. 2.0, NASA RP- 1190. 



This object is barely visible in the IRAC 4.5 /xm image, and 
it is not seen in the 2MASS images of the field. The IR-radio 
spectrum and the lack of an extended optical or NIR counter- 
part suggest that the IR emission is thermal emission of dust 
with T ~ 30-50 K, possibly mixed with the molecular/atomic 
cloud(s) that are heated and partly ionized by radiation from 
embedded young stars. Its luminosity can be crudely esti- 
mated as La ~ 3 x 10 37 c/| ergs s" 1 . 

The nature of the other two bright radio sources within the 
brightest part of HESS J1825 (B and C in Fig. 3f), centered 
at a » 18 h 25 m 57 s , 5 w -13°47'56" and a » 18 h 25 m 38 s , S w 
-13°36'37", respectively, is unknown. However, as they are 
not seen in IR, their bolometric luminosity is much lower than 
that of Source A. 

3. PROPER MOTION OF PSR B1823-13 

We measured the proper motion using two unrelated 
archival VLA observations. Both observations were made 
in the 20 cm band in the A-configuration providing about 
1.5" x 2.1" resolution. The reduction technique is identi- 
cal to that used by Brisken et al. (2006) in the determina- 
tion of the proper motion of PSR B 1800-21 using 1.4 GHz 
data. All synthesis calibration and imaging were performed 
using the AIPS software package. Project AF118 was ob- 
served on MJD 46573 in right circular polarization in spec- 
tral line mode providing 15 channels spanning 50 MHz band- 
width at 1385 MHz. The upper 4 channels were discarded 
due to severe interference. Project AC761 was observed on 
MJD 53265 in dual circular polarization in spectral line mode 
providing 15 channels across 25 MHz at 1665 MHz. After 
standard calibration and editing, AIPS task UVFIX was used 
to correct the baseline vectors of the pulsar field for the half- 
integration time time-stamp error common to all VLA obser- 
vations as well as including relativistic aberration correction. 
UVFIX was also used to convert the B1950 coordinates used 
in AF1 18 to the J2000.0 coordinates. 

The temporal and spectral resolution of both data sets al- 
lowed imaging of the entire 25' primary beam. In addition to 
the pulsar, eleven other sources were detected in both epochs, 
five of which were bright and compact enough to be used 
as reference sources. Accurate positions for the pulsar and 
these five objects were measured for each epoch. The motion 
of the pulsar was then measured in a frame defined by these 
sources using a fitting prodedure described in McGary et al. 
(2001), yielding a proper motion /j, a = 22.7 ±2.2 mas yr" 1 , 
fig =-5.3 ±2.7 mas yr" 1 (68% confidence intervals). The 
global linear fit used to determine the proper motion and a 
possible translation and rotation of the coorinates between 
the two epochs yielded a reduced \ 2 of 1.3, indicating that 
the formal uncertainties may be underestimates at the 15% 
level. Assuming a distance of 4 ± 1 kpc, this proper mo- 
tion implies a local standard of rest (LSR) proper motion of 
fi a = 23 .0 ± 2.5 mas yr" 1 , = -3 .9 ± 3 . 3 mas yr" 1 , which cor- 
responds to the velocity v± = (443 ±48)d4 km s" 1 at a position 
angle of 99.°6±8.°1 (east of north). At this proper motion, the 
pulsar has traveled about 8.2' in the plane of the sky during its 
characteristic age of 21 kyr. The projected birthplace of the 
pulsar (hence the SN location) is marked in Figure 3b-f. 

4. DISCUSSION 

4.1. Pular Wind Nebula 

From their analysis of the XMM-Newton data, G03 have 
concluded that the B1823 PWN (which they designated as 
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G18. 0-0.7) is comprised of two components ("core" and "dif- 
fuse") of different brightness and size and perhaps different 
spectra. Our Chandra ACIS observation has confirmed the 
presence of the two components, which we designate as the 
"compact PWN" and "extended PWN", and allowed us to in- 
vestigate the properties of the compact PWN in much more 
detail. The better resolved images and the measurement of 
the pulsar's proper motion not only allow a more definitive 
interpretation of the compact PWN properties, but they also 
help better constrain the models the puzzling extended PWN 
and its connection to HESS J1825. 

4.1.1. Compact PWN 

The east-west elongation of the compact PWN, approxi- 
mately coinciding with the direction of the pulsar's proper 
motion (see Fig. 2), strongly suggests that the compact PWN 
is shaped by the oncoming flow of the ambient matter in 
the pulsar's reference frame. The transverse velocity of the 
pulsar, vj_ w 440c/4 km s" 1 , is high enough to assume that 
the pulsar's motion is supersonic, i.e., M. = v/c s > 1, or, 
equivalently, T < 1.4 x 10 7 p,(d 4 / sin i) 2 K, where M is the 
Mach number, v = vj_/ sin; is the pulsar's speed, i is the an- 
gle between the velocity vector and the line of sight, c s = 

(5kT /3/xm H ) 1//2 = lllfT l ^ 2 T^' 2 km s" 1 is the speed of sound, 
p is the molecular weight, and T = lO^K is the temperature 
of the ambient medium. When a pulsar moves supersonically, 
the pulsar wind (PW) termination shock (TS) acquires a bow- 
like shape ahead of the pulsar, with the apex at a distance 
Rh = (£/h/47rcp ram ) 1//2 , determined by the balance between 
the PW pressure and the ram pressure, p mm = pv 2 , where p is 
the density of the ambient medium, and fa takes into account 
anisotropy of the pulsar wind. For B1823, these equations 
yield p mm = 3.3 x \Q~ 9 n{di,/ sin;') 2 ergs cm" 3 and Rh = 4.8 x 

10 16 n _1 / 2 /^ 2 (c/4/sin0 _1 cm, which translates into the pro- 
jected angular distance Hh « 0.%" fl{ 2 n~ l l 2 (dn/ sin/)" 2 , where 
n = p/mn. 

In an idealized picture, the shocked PW, responsible for 
the PWN X-ray emission, is confined between the TS and 
the contact discontinuity (CD) surface, which separates the 
shocked PW from the shocked ambient medium. For A4 3> 1 
and a nearly isotropic preshock PW with a low magnetization 
parameter a, the distance between the pulsar and the head of 
the CD surface is w 1.3/?/, (Bucciantini et al. 2005, hereafter 
B05). In a more realistic picture, one should not expect a 
smooth CD surface because it is unstable. As a result, the 
shocked PW can penetrate into the shocked ambient medium, 
and vice versa. 

The images of the compact PWN (Figs. 1 and 2) do not 
show a sharp PWN boundary ahead of the pulsar that could 
be unambiguosly identified as a CD surface head. Instead, 
we see a two-component structure, with the bright inner com- 
ponent up to 3" east of the pulsar and some fuzzy emission 
discernible up to 7". Although the quality of the image does 
not allow one to make a definitive conclusion, we can as- 
sume that the fuzzy emission represents the shocked PW that 
has penetrated into the shocked ambient medium because of 
CD instabilties, while the inner component is emitted by the 
shocked PW immediately beyond the TS, so that IZi, ~ 2" -3" 
seems to be a reasonable estimate for the angular distance to 
the TS head. Scaling TZf, to 3", we can estimate the ambient 
number density, n ~ 0.07/h(c/ 4 / sin0~ 4 (7^/ ! /3")~ 2 cm" 3 , and 
the ram pressure, p mm ~ 2.3 x 10" 1 %(^4/sin/)- 2 (7e ft /3")" 2 



ergs cm -3 . Such a ram pressure is typical for Vela-like pul- 
sars (Kargaltsev et al. 2007), while such a number density can 
be found in a warm phase of the ISM or within an evolved 
SNR. As no SNR emission is seen in the vicinity of B1823, 
we cannot measure the temperature and pressure of the am- 
bient medium, but it seems reasonable to assume T ~ 10 5 K, 
which could explain the lack of X-ray SNR emission. Scaling 
the temperature to this value, we obtain p. dm \, = nkT //j~lx 

\Q- n fn{d,/^m)-\n h /yr 2 ^ X T 5 < Pmm , c s ^ 40p-^ 2 T^ 2 

km s~\ and M ~ 10( 1 i/siniV /2 7r 1/2 . 

The existing MHD models of bow-shock PWNe associated 
with supersonically moving pulsars assume a nearly isotropic 
preshock PW. However, the high-resolution Chandra images 
of the whole compact PWN can hardly be reconciled with this 
assumption. According to B05, the TS in an isotropic PW of 
a supersonically moving pulsar acquires a bullet-like shape, 
with the bullet's cylindrical radius r TS ~ R h and the distance 
of its back surface from the pulsar R/, ~ 6Rh- The shocked 
PW outside the TS is confined inside the CD surface, which 
has a cylindrical shape behind the TS bullet, with a cylindrical 
radius rco ~ 4/?/,. The collimated PW flows with subrelativis- 
tic velocities: 0.1c-0.2c in the inner channel, r < r T s, and up 
to 0.8c-0.9c in the outer channel, rxs < r < rcD (see Figs. 
1-3 in B05). The images of the compact PWN are not quite 
consistent with this simple picture. First of all, the bright in- 
ner component, which could be naturally associated with the 
PW immediately outside the TS, is inclined to the direction 
of the proper motion by an angle <~ 50°, while it should be 
coaligned with the proper motion if the wind is isotropic. Sec- 
ond, the outer PWN, which might represent the shocked PW 
confined by the CD surface, not only has an irregular bound- 
ary (which could be explained by instabilities), but it is also 
rather nonuniform. (Particularly puzzling is the blob m 8" 
west of the pulsar; see Figs. 1 and 2, and §2.2). Finally, we do 
not see an extended tail behind the TS aligned with the proper 
motion direction. These differences between the model and 
the observed image suggest that the preshock PW is, in fact, 
substantially anisotropic. To understand the PW geometry 
and interpret the PWN morphology more or less unambigu- 
ously, deeper high-resolution X-ray observations are needed, 
as well as modeling of intrinsically anisotropic, magnetized 
PWs from fast-moving pulsars. Here we can only speculate 
that the inner PWN likely represents a shocked equatorial out- 
flow (similar to the famous Crab torus; Weisskopf et al. 2000), 
seen almost edge-on. If this is the case, then the spin axis of 
the B1823 pulsar is strongly misaligned (by ~ 40°-50° in the 
plane of the sky) with the direction of its motion, contrary to 
many other pulsars. The absense of jets along the spin axis in 
the observed images can be explained by their faintness (e.g., 
the jets of the Vela PWN would be unobservable at d ~ 4 kpc 
with such an exposure - see Fig. 8 in Kargaltsev et al. 2007), 
or possibly the jets are deflected (or even crushed) by the on- 
coming flow of the ambient medium. The blob and the appar- 
ent wiggles in the west part of the outer component could be 
ascribed to instabilities in the tail formed by the supersonic 
motion of the pulsar. 

The luminosity, L comp ~ 3 x 10 32 c/ 2 ergs s" 1 in the 0.5-8 
keV band 8 , the X-ray efficiency, ?7 com p = L comp /E w 1.2 x 

8 We should note that our estimate for the compact PWN luminosity is a 
factor of 3 lower than that obtained by G03 for the "core" component. The 
discrepancy is possibly caused by the much larger area from wich G03 ex- 
tracted the spectrum (1018 vs. 225 arcsec 2 ), which included some contribu- 
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10~ 4 <f|, and the spectral slope, r comp w 1.3 (at the assumed 
«h,22 = 10) of the compact PWN are not unusual for Vela- 
like PWNe (e.g., Kargaltsev et al. 2007; Kargaltsev & Pavlov 
2007). Also, similar to other Vela-like PWNe, the PWN lu- 
minosity exceeds by a factor of a few the pulsar luminosity in 
the same X-ray energy range, L comp ~ 4L P sr. 

4.1.2. Extended PWN 

The morphology of the extended PWN is even more puz- 
zling than that of the compact one. In the conventional picture 
of a bow-shock PWN, we would expect to see an extended 
component behind the compact PWN of the moving pulsar, 
which would look like a "tail". In B1823, however, the PWN 
extends southward from the pulsar, almost perpendicular to 
the proper motion direction, similar to the extended diffuse 
PWN of the Vela pulsar (Pavlov et al. 2003). This means that 
something is deflecting (or has deflected) the shocked PW in 
this direction. One could imagine a local north wind in the tur- 
bulent ambient medium, but such a wind would have to blow 
with a velocity of several hundred km s" 1 (comparable to the 
pulsar velocity) to displace the PWN bubble. If such a wind 
were blowing now, we would see the compact PWN elongated 
southward from the pulsar, and, more importantly, it would be 
not a usual wind but a shock, since the required velocity al- 
most certainly exceeds the sound speed. Thus, it seems we 
have to assume that the PWN bubble was crushed/displaced 
by a passing shock some time ago. A natural candidate for 
such a shock is the reverse shock from the surrounding SNR 
that has propagated back inward and collided with the PWN, 
displacing the PWN bubble from the pulsar, as has been sug- 
gested by Blondin et al. (2001) for the strongly dispaced Vela 
X component of the Vela PWN. Such an explanation for the 
one-sided extended PWN of B1823 has been proposed by 
G03, and we currently do not see an alternative explanation. If 
this hypothesis is correct, we may expect the extended PWN 
to consist of two components: a strongly displaced one (simi- 
lar to Vela X), which represents a relic population of relativis- 
tic electrons that have lost a substantial fraction of their en- 
ergy to synchrotron, inverse Comption, and adiabatic cooling, 
and a younger one, which represents electrons that entered the 
PWN after the inverse shock passage. 

Although the crushed/dragged PWN hypothesis looks at- 
tractive and provides a plausible explanation for the extended 
X-ray PWN, we should mention that it does not so readily 
expalin the huge, about 70 pc at d = 4 kpc, size of the TeV 
PWN. Traveling 70 pc in time t requires the average velocity 
v w 6,900(f/10kyr) _1 km s" 1 , too high for a reverse shock. 
Even if we assume t = 21 kyr (i.e., equal to the pulsar's spin- 
down age), we still need v w 3,300 km s" 1 . Therefore, it is 
hard to imagine how the reverse shock could have dragged 
the relic PWN so far, even if it passed the pulsar very long 
time ago. Perhaps, one should explore another possibility, that 
during a substantial fraction of the pulsar's lifetime the PWN 
electrons have been transported southward from the pulsar 
along an anisotropic large-scale magnetic field. Alternatively, 
one might assume a strongly anisotropic supply of electrons 
into the PWN, e.g., through a one-sided jet (which, however, 
is not seen in the images available). At this point, however, 
these possibilities remain highly speculative. 

tion from the extended ("diffuse") PWN, the pulsar, and field sources, such as 
Source X in our Figures 1 and 2. Also, G03 extracted the background spec- 
trum from a region outside the extended PWN, while we used an annulus 
filled with the extended PWN emission (see §2.3). 



As has been noticed by G03, the extended PWN spectrum 
is softer than the compact PWN spectrum, which is also con- 
firmed by our results (see Fig. 5). The difference of the pho- 
ton indices, r ext -r comp w 0.5, supports the hypothesis that 
the extended X-ray PWN electrons have lost their energy to 
synchrotron cooling after traveling a distance of a few par- 
sees. The slope r ext ~ 1 .9 ± 0. 1 we found for the X-ray spec- 
trum is very close to the slope of the VHE 7-ray spectrum 
within w 6' from the pulsar, T 7 = 1.8-2.0 (A06), which con- 
firms that both the X-ray and 7-ray emission are due to the 
same electron population with a power-law energy spectrum 
dN e /dE e oc E~ p , with p = 2r ext - 1 « 2.5-3. Interestingly, the 
TeV spectrum at larger distances from the pulsar, where we 
do not see the X-ray PWN, becomes even softer (T 7 grows up 
to w 2.6 at 60' from the pulsar). 

Because we can see only part of the extended PWN in our 
ACIS image (see §2.2), we can measure only a fraction of 
its X-ray luminosity (L ext ~lx 10 33 d 2 ergs s" 1 in the 3.4 
arcmin 2 region shown in Fig. 1, bottom right). The total ex- 
tended PWN luminosity should be a factor of a few higher. 
For instance, G03 found L ext ~3x I0 33 dj ergs s" 1 in a 33.4 
arcmin 2 area, which corresponds to the extended PWN effi- 
ciency f7 ext ~ 10" 3 c/|, a factor of ~ 10 higher than that of the 
compact PWN. 

4.2. Pulsar 

Although we have resolved the pulsar from the PWN, the 
small number of pulsar counts detected and the strong and not 
well determined background from the inner component of the 
compact PWN preclude a precise spectral analysis. We can 
only estimate the pulsar luminosity, L P sr ~ 1 x 10 32 dj ergs s" 1 
in the 0.5-8 keV band, which is a small fraction, ~3x 10" 5 , 
of the pulsar's spindown power E. The slope of the pulsar 
spectrum, Tpsr ~ 2-3, is strongly correlated with the ISM 
column density (see Fig. 7). It is likely that the spectrum con- 
sists of two components, magnetospheric and thermal, but the 
photon statistics does not allow one to separate these compo- 
nents unambiguously. For «h.22 = 1 .3 (at the higher end of the 
plausible «h range) and projected emitting areas in the range 
of (20-500)4 km2 ' we found kT ~ 90-150 eV and L bol ~ (3- 
20) x 10 32 d 2 ergs s" 1 for a possible BB component (higher 
luminosities correspond to lower temperates and larger emit- 
ting areas). These temperatures and luminosities are generally 
similar to those of other Vela-like pulsars (e.g., kT w 130 eV, 
A w 20(t//300pc) 2 km 2 , L bol - 2 x 10 32 (t//300pc) 2 ergs s" 1 
for the Vela pulsar; Pavlov et al. 2001). However, because of 
the large uncertainties and poorly constrained n H , we cannot 
be sure even in the detection of the thermal component, not 
to mention measuring its parameters. Therefore, we should 
treat the estimated bolometric luminosities of B 1823 as upper 
limits rather than actually measured values. The most conser- 
vative of the above limits, Lbol ~ 2 x 10 33 d 2 ergs s" 1 , corre- 
sponding to A = 500<f| km 2 and T « 100 eV, is close to the 
bolometric luminosities predicted by the so-called "standard" 
NS cooling model for the NS age of 20 kyr (see, e.g., Fig. 2 
in Tsuruta 1998). Overall, we can conclude that B1823 looks 
like a typical Vela-like pulsar, but a deeper high-resolution 
observation is needed to infer its properties in more detail. 

4.3. The lack of a host SNR? 

We see no emission in our ACIS images that could be as- 
cribed to an SNR associated with B1823, and no such emis- 
sion has been reported from the previous X-ray (F96; G03) or 
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radio (Braun et al. 1989; Brogan et al. 2006) observations. Us- 
ing our proper motion measurement, we extrapolated the pul- 
sar trajectory back in time to check whether or not it crosses 
a known SNR. The nearest SNR to the pulsar trajectory is 
the 8' x 8' shell SNR G18. 16-0.16 (Brogan et al. 2006) at a 
distance of 35' from the current pulsar position. Since the ex- 
trapolated trajectory misses the SNR center by 19' ±4', we 
conclude that G 1 8 . 1 6-0. 1 6 is not related to B 1 823 . 

The failure to detect the host SNR in X-rays can be at- 
tributed to the strong ISM absorption of the soft thermal emis- 
sion from the SNR interiors (see Kargaltsev et al. 2007, where 
a similar situation is discussed regarding the Vela-like PSR 
B 1800-21, which also lacks an associated SNR). Therefore, it 
is possible that B 1823 is still within the interiors of an evolved 
SNR, which is undetectable in X-ray or radio, but whose pres- 
sure can contribute to the PWN confinement. 

4.4. Connection between B1823 and HESS J 1825 

The important discovery of the VHE 7-ray spectral soften- 
ing with increasing distance from B1823 (A06) has provided 
further evidence that this pulsar/PWN is the main source of 
the ultrarelativistic electrons powering the TeV plerion. Qual- 
itative models for this VHE 7-ray source have been discussed 
in a number of works (A06; de Jager et al. 2005; de Jager 
2007). Generally, they interpret the 7-ray emission as pro- 
duced by the Compton upscattering of the CMB photons by 
electrons with energies E e ~ 25(£ 7 /TTeV) 1/2 TeV - 10-100 
TeV [electron Lorentz factor j e ~ (2-20) x 10 7 ], supplied by 
the B1823 pulsar. The fact that the pulsar is located at the 
northern boundary of the TeV plerion is explained by the ef- 
fect of the SNR reverse shock that has dragged the relic PWN 
southward. (An alternative explanation, that the TeV plerion 
represents electrons left behind the pulsar moving northward, 
is invalidated by our measurement of the proper motion.) The 
softening of the TeV spectrum with the distance from the 
pulsar can be ascribed to the radiative losses (mostly to the 
IC scattering losses) or some energy-dependent "convection", 
or variation of the particle injection spectrum with the pul- 
sar's age. As A06 note, even at very low ambient magnetic 
fields, B < 3 jUG, when the IC energy losses dominate the syn- 
chrotron losses, the 7-ray efficiency, 7/ 7 = L 7 /E w 0. 1 , is unre- 
alistically high if we use the current value of E. This supports 
the hypothesis that the VHE 7-ray emission is produced by 
relic electrons emitted by the pulsar when it was very young 
and much more powerful. 

Modelling of the TeV plerion is beyond the scope of this 
work. We only note some open problems with the current in- 
terpretation of this remarkable object. First, as we have men- 
tioned in §4.1, it is difficult to reconcile the very large size of 
the TeV plerion with the hypothesis that the relic PWN was 
dragged by the SNR reverse shock to its current position. In 
principle, other sources of ultrarelativistic electrons within the 
1° size field could contribute to the observed TeV emission 
and increase the apparent size of the TeV source. Our anal- 
ysis of several X-ray sources in the brightest part of the TeV 
plerion has not shown any obvious candidates (see §2.5), but 
one cannot rule out the possibility that there are such sources 
outside the small ACIS field-of-view. 

Another possible problem with the current interpretation of 
the TeV source is the spatial distribution of its brightness and 
spectrum. The brightest part of the TeV plerion is offset by 
~ 10' (~ 12 pc in the plane of the sky) south-southwest of the 
pulsar (see Fig. 1 of A06), while the plerion extends ~ 50' fur- 



ther in the same direcion (almost perpendicular to the direc- 
tion of the pulsar's proper motion). One could speculate that 
the relic PWN was displaced by ~ 10' some time ago [which 
would require a reverse shock speed of ~ 1000(f/10kyr) _1 
km s" 1 ], and its relativistic electrons have been diffusing from 
that position. It is however, unclear as to why they diffused 
not radially but preferentially in the south-southwest direc- 
tion. If the passing shock not just displaced the relic PWN to 
a certain distance but smeared it all over the whole extent of 
the currently seen TeV plerion, then a very high shock speed 
would be required (see §4.1). Furthermore, the spectral map 
of HESS J 1825 (see Figs. 5 and 6 of A06) shows that the VHE 
7-ray spectrum softens radially not from the maximum of the 
brightness distribution but from the pulsar position. This hints 
that at least some contribution to the TeV photons comes from 
the same young electrons that are responsible for the extended 
X-ray PWN of a much smaller size. This conjecture is fur- 
ther supported by the fact that the spectral slopes of the X-ray 
(synchrotron) and 7-ray (IC) spectra are about the same at a 
few arcminute distance from the pulsar, which suggests that 
they are generated by the same population of electrons. On 
the other hand, if the TeV 7-ray emission is due to the IC 
scattering of the CMB photons, and the relativistic electrons 
are mostly supplied by the pulsar, one should expect the TeV 
brightness to peak at the pulsar position, not 10' from the pul- 
sar. This inconsistency might be alleviated assuming that in 
the immediate vicinity of the pulsar, where the electrons are 
more energetic, they upscatter the CMB photons in the Klein- 
Nishina regime, which results in a lower brightness of the IC 
scattered radiation, but it remains to be seen from detailed 
models whether or not this assumption is realistic. Overall, it 
seems plausible at this point that, at least in the extended X- 
ray PWN region, the TeV emission can be produced by two 
electron populations of different ages. To infer the properties 
of these populations, deep observations with XMM-Newton 
and GLAST of a region from which both X-ray and VHE 7- 
ray emission is seen would be very helpful. 

Another issue related to the nature of the TeV plerion and its 
connection to the B 1 823 pulsar/PWN is whether the TeV pho- 
tons are produced by upscattering of the CMB or some other 
radiation. Close to the Galactic plane, where HESS J 1825 is 
situated, a factor of a few higher radiation energy density can 
be provided by IR emission from interstellar dust and Galac- 
tic starlight. For instance, at the galactocentric distance of 
4 kpc, the models of interstellar radiation by Strong et al. 
(2000) give the energy densities of 0.6 and 2.7 eV cm" 3 for 
these two components, respectively, versus 0.26 eV cm -3 for 
the CMB. Moreover, if there are very luminous sources of 
radiation (such as molecular/dust complexes in star-forming 
regions) in the vicinity of the source of ultrarelativistic elec- 
trons, the high radiation energy density around these sources 
can strongly enhance the TeV luminosity (see, e.g., Kargalt- 
sev & Pavlov 2007). In this regard, the bright radio-IR Source 
A (see Fig. 3), whose spectrum peaks at A > 100 /im, might 
be of some interest. With the bolometric flux Fa ~ 2 x 10~ 8 
ergs cm" 2 s -1 , its radiation energy density exceeds the am- 
bient radiation energy density U dm b at angular distances 9 < 
(F A /cU- dbm ) 1/2 ~ 2.2'(t/ amb /leVcm- 3 r 1/2 . If the distance to 
Source A is about the same as to B1823, and the TeV emis- 
sion is produced by the IC scattering of IR photons with en- 
ergies around 0.03 eV (which would require electron energies 
E e <~ 1-20 TeV), then Source A could create a few arcminute 
size brightening. Moreover, if there are hot and massive stars 



Chandra Observation of PSR B 1823-1 3 



11 



hidden within this dust/molecular complex of ~ Ad$ pc size, 
their shocked winds could be additional sources of TeV elec- 
trons, further enhancing the VHE 7-ray emission from this 
region. (VHE 7-rays apparently produced by stellar winds in 
the open stellar cluster Westerlund 2 have recently been de- 
tected by Aharonian et al. 2007.) Although this hypothesis 
looks rather speculative at this point, and Source A can not be 
responsible for the bulk of the detected TeV radiation, ded- 
icated observations of this region in the IR and radio would 
help evaluate the Source A contribution to HESS J1825. 

5. CONCLUSIONS 

The high-resolution Chandra ACIS observation has allowed 
us to resolve the B1823 pulsar from its PWN, image the PWN 
structure, and measure the spectra of the pulsar and the PWN 
components. We found that the pulsar's X-ray luminosity, 
Lpsr ~ 0.8 x 10 32 ergs s -1 in the 0.5-8 keV band, and spec- 
trum are similar to those of other Vela-like pulsars, but the 
small number of pulsar counts and the contamination from the 
inner PWN component do not allow an unambiguous separa- 
tion of the thermal and nonthermal components of the pulsar 
radiation. 

We have measured the pulsar's proper motion and found its 
transverse velocity, v_i_ = 443 ±48 km s" 1 at d = 4 kpc, close to 
an average velocity of young radio pulsars. We found no in- 
dication of an SNR, in neither X-rays nor radio, between the 
current position of the pulsar and its birthplace, ~ 8' west- 
ward. 

We have confirmed the presence of the compact and ex- 
tended PWNe of different shapes, surface brightnesses, and 
spectra. The compact PWN, of ~ 0.5 x 0.2 pc 2 size in the 
plane of the sky, is elongated approximately along the proper 
motion direcion, with the pulsar close to its leading edge. This 
suggests that the postshock pulsar wind in the compact PWN, 
responsible for its X-ray emission, is confined by the ram 
pressure due to the supersonic motion of the pulsar. The X- 
ray spectrum and luminosity, L comp ~ 3 x 10 32 ergs s" 1 , of the 
compact PWN are similar to those of other Vela-like PWNe. 

The X-ray surface brightness distribution within the com- 
pact PWN is very nonuniform. In particular, we see a much 
brighter inner PWN component, ~ 0.14 x 0.06 pc 2 , around 
the pulsar, inclined by ~ 50° with respect to the direction of 
the proper motion. Such a structure indicates that the pulsar 
wind is intrincically anisotropic. If the inner PWN compo- 
nent represents the shocked equatorial outflow from the pul- 
sar magnetosphere, the pulsar's spin axis is strongly inclined 
to the proper motion direction, in contrast to many other pul- 
sars. 

The extended PWN is mostly seen south-southwest of the 
pulsar (up to 2.4' from the pulsar, in the ACIS image), almost 
perpendicular to the proper motion direction. Its softer spec- 
trum indicates that its X-ray emission is generated by older 
electrons that have lost part of their energy because of the ra- 
diative cooling. The southward displacement of the extended 
X-ray PWN with respect to the pulsar might be caused by 
a strong pressure gradient in the ambient medium and/or by 
the reverse SNR shock, which has swept the PWN material 
southward. 

The extended X-ray PWN is apparently connected to the 
VHE 7-ray PWN that extends in the same direction but to a 
much larger distance from the pulsar, up to 70 pc in the plane 
of the sky. The 7-ray emission is likely due to the IC scatter- 
ing of background radiation off the relativistic electrons with 
energies up to 100 TeV supplied by the B1823 pulsar. How- 



ever, the very high luminosity, L 7 ~ 3 x 10 35 ergs s -1 , can only 
be explained if the bulk of 7-radiation is generated by relic 
electrons, produced long ago, when the pulsar was much more 
powerful. The relic electrons could be displaced from the pul- 
sar by a reverse SNR shock, but it remains unclear as to why 
the 7-ray PWN is so large. We have not found X-ray sources 
in the brightest part of the 7-ray PWN which could signif- 
icantly contribute to the 7-radiation. Although the IR-radio 
images of the field show a star-forming region in the bright- 
est part of the 7-ray PWN, which could enhance the 7-ray 
emission by supplying additional photons for the IC scatter- 
ing, and perhaps even additional relativistic electrons, it can 
be responsible only for a small fraction of the observed VHE 
7-radiation. To better understand the multiwavelength nature 
of the complex B1823 PWN, deep X-ray and IR observations 
of the field would be particularly useful. 
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TABLE 1 
PL FITS TO THE PWN SPECTRUM 

Region A/"" F F b F unabs c X 2 /do~ 

Compact PWN 2. 14 ±0.25 1.33^ ^ 1.1 ±0.1 1.7±0.1 0.93/18 
ExtendedPWN 9.0±1.0 1.87±0.14 2.5±0.2 4.5;°; 3 1.28/13 

NOTE. — The fits are for fixed «h,22 = »h/(10 22 cm" 2 ) = 1.0. The uncertainties are given at 68% confidence level for a single interesting parameter. 
"Spectral flux in units of 10" 5 photons cm " 2 s" 1 keV" 1 at 1 keV. 
b Observed energy flux in units of 10" 13 ergs cm" 2 s _1 , in the 0.7-7 keV band. 
c Unabsorbed energy flux in units of 10" 13 ergs cm" 2 s in the 0.5-8 keV band. 



TABLE 2 
Fits to the pulsar spectrum 



Model 


"H.22 


A/" a or A b 


T or kT c 


C 


L x or L bo i d 


PL 
PL 
PL 

PL+BB(PL) 


0.7 
1.0 
1.3 
1.3 


19.6 + |5 

7 5+ 3 - 6 
'•-"-l .9 


1 98+ - 40 
9 4^+0.53 

1.9±0.6 


639 
646 
657 
650 


0.6 
0.8 
1.1 
0.7 


PL+BB(BB) 


1.3 


500 


97+4 


650 


18 


PL+BB(PL) 


1.3 


5.0+5; 1 


1.7 ±0.7 


651 


0.6 


PL+BB(BB) 


1.3 


20 


139^ 


651 


3.0 



NOTE. — The fits are for fixed «h.22 = "h/'O 22 cm 2 . The uncertainties are given at 68% confidence level for a single interesting parameter. In each case the 
fits were done using the C-statistics (Cash 1979) and the unbinned source and background spectra with the total of 1052 channels. 
"Spectral flux in units of 10" 6 photons cm " 2 s" 1 keV" 1 at 1 keV. 

b Projected area of the emitting region for the BB model, in units of km 2 , fixed in the fits. 
C BB temperature in eV. 

''Unabsorbed PL luminosity in the 0.5-8 keV band or bolometric BB luminosity, in units of 10 32 ergs s _1 . 



TABLE 3 

Positions and spectral properties of the field X-ray sources 

Source" R.A. Peel. 5r b n H .22 A/~ c F F ob s J jWj? 

X 18:26:14.43 -13:34:48.0 0.4" 2.9+ 2 |? 0.10+jj^ 0.2±0.5 2.9±0.5 3.7 

si 18:26:14.06 -13:40:35.8 0.5" 4.3 +i: ' 2.2+\ J 2.1 ±0.5 3.2±0.4 11.9 

s2 18:26:17.18 -13:41:12.7 0.5" 0.9^ 1.8^ 3.1^ 3 1.5 ±0.3 7.0 

s3 18:26:17.99 -13:42:27.5 0.5" 0.5+{Jj 0.5±0.2 1.9±0.4 1.7±0.4 2.7 

s4 18:26:20.88 -13:44:25.7 0.8" 1.2*%* 3.5^ 3.1 ±0.3 2.8±0.3 14.2 

s5 18:26:24.76 -13:37:06.2 0.5" 3.0+f« 0.8+{J : « 1.8±0.5 1.9±0.3 4.6 

il 18:25:40.55 -13:40:11.2 0.6" l.l+j; 2 , 1.3+gj 1.9±0.3 3.2±0.6 6.4 

NOTE. — The uncertainties of the fitting parameters are given at 68% confidence level for a single interesting parameter. 
"Source designation as shown in Figs. 1 and 8. 
b Radius of the position error circle at the 68% confidence level. 
c Spectral flux in units of 10~ 5 photons cm " 2 s"' keV" 1 at 1 keV. 
d Observed flux in the 1-8 keV band, in units of 10~ 14 ergs cm" 2 s~' . 

e Best-fit value of the unabsorbed PL flux in the 0.5-8 keV band, in units of 10"' 4 ergs cm" 2 s _1 . 



